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INTRODUCTION 


When  a step-function  shock  wave  is  introduced  into  a medium, 
with  the  result  that  a multi-step  function  disturbance  propagates  with- 
in the  medium,  shock  wave  physicists  refer  to  that  medium  as  unstable. 

A medium  can  also  be  unstable  with  respect  to  fluctuations  (e.g., 
volume  fluctuations)  in  the  absence  of  an  externally  induced  dynamic 
strain.  This  type  of  instability  is  called  STATE  INSTABILITY.  In  order 
for  a multi-step  shock  structure  to  be  observed,  each  step  state  must 
be  state  stable.  This  report  will  treat  aspects  of  both  shock  and  state 
stability  and  instability. 

Shock  wave  stability  studies  usually  involve  transitions  between 
known  stable  states.  This  study  is  interested  in  the  possibility  of 
shock  transitions  to  new  and  relatively  exotic  media  and  is  therefore 
concerned  with  b6th  types  of  stability.  In  the  process  of  this  dual 
approach,  it  is  shown  that  the  requirements  for  shock  instability  and 
state  stability  (i.e.,  the  longtime  existence  of  a state  achieved  via  a 
hypothetical  shock  transition)  can  be  mutually  exclusive  in  terms  of 
materials  properties. 

One  of  the  most  interesting  aspects  of  shock  instability  theory 
involves  the  possibility  of  shock  amplitude  amplification.  Such 
amplification  can  occur  (ref.  1)  when  small  amplitude  acoustic 
disturbances,  originating  behind  a propagating  shock  front,  catch 
up  to  and  reflect  from  the  impedance  discontinuity  of  that  front.  If 
materials  properties  are  correctly  chosen  (corresponding  to  shock 
instability  criteria) , the  acoustic  reflection  amplitude  can  be  larger 
than  the  amplitude  of  the  prereflection  acoustic  signal . Such  a 
phenomenon  represents  amplification  of  the  acoustic  signal  and  results 
in  amplification  of  the  amplitude  of  the  shock  front. 

Acoustic  signals  originating  behind  the  shock  front  can  be  obtained 
via  a member  of  mechanisms: 

1 . Local  exothermic  chemical  reaction  behind  the  shock  front 
can  send  out  differential  compressive  disturbances  via  a thermal  ex- 
pansion coefficient. 
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2.  Grading  (i.e.  varying)  the  impedance  properties  of  an  im- 
pacting projectile  will  produce  a ramp  disturbance  which  can  be  approxi- 
mated by  an  acoustic  disturbance  riding  on  a step  shock. 

3.  A second  order  (isovolumetric  without  latent  heat)  phase 
transition,  which  results  in  an  elastically  harder  material,  will  lead  to, 
for  a given  strain,  a higher  stress.  Thus,  if  such  transitions  occur 
behind  the  shock  front,  the  resulting  higher  stress  regions  will  propagate 
as  acoustic  disturbances. 

In  order  for  the  amplification  to  result  in  avalanching,  for  one  incoming 
acoustic  disturbance,  the  signal  which  is  reflected  from  the  shock  front 
must  lose  less  amplitude  at  the  impacted  surface  than  was  gained  upon  re- 
flection from  the  shock  front.  This  is  essentially  a statement  about  the 
relative  impedance  of  the  driver. 

Regardless  of  whether  one  has  single-pass  amplification  or  amplification 
with  some  avalanching,  the  amplification  process  represents  the  possibility 
of  using  materials  properties  to  pulse-shape  a shock  front.  This  is  es- 
pecially exciting  when  the  shocked  medium  is  inert. 

If,  in  an  inert  material , a shock  of  amplitude  P with  pulse  width  t is 
introduced  and  pulse  shaped  to  amplitude  2P  and  pulse  width  t/2,  then  the 
pulse-shaped  shock  disturbance  would  be  more  efficient  than  the  original 
shock  in  initiating  detonation  in  explosives.  A commonly  used  criteria 
(ref.  2)  for  detonability  is  P2t  = constant  so  that  the  pulse-shaped  shock 
would,  in  a sense,  be  more  efficient  by  a factor  of  two.  Other  applications 
exist. 

SIMPLE  SHOCK  WAVE  INSTABILITIES 

Exceeding  the  Hugoniot  elastic  stress  limit  (ref.  3)  or  exceeding  the 
stress  necessary  for  a polymorphic  (structural)  phase  transition  (ref.  3 
and  4)  (e.g. , a toe  iron  at  approximately  130  kbars)  represent  the  two 
most  common  experimental  methods  of  generating  the  multi-step  shock 
structure  characteristic  of  a medium  exhibiting  a shock  instability.  In 
such  simple  cases  the  P-V  and  P-u  curves  appear  as  shown  in  figures  1 

and  2. 


Figure  1.  P-V  Hugoniot  segments  for  elastic  limit  or  polymorphic 
transition  occurring  at  VA. 
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Figure  2.  P-u  Hugoniot  segments  for  elastic  limit  or  polymorphic 
transition  occurring  at  VA. 
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The  important  point  of  figures  1 and  2 is  the  sign  of  — | 


H 


. dP  . 

andsrJ  h 


(the  subscript  H denoting  Hugoniot*)  for  each  of  the  various  segments.  For 
these  simple  shock  instabilities  it  will  be  shown  that 


dV1  H 


< 0 


. dP 
and  — 
du 


„ > 0 ALWAYS  hold  . 
H 


For  a steady  state  shock  state  B propagating  into  the  preshocked  state  A, 
the  conservation  of  mass  and  momentum  can  be  expressed  in  the  jump  con- 
ditions below: 


u - uA  = PA  (U  - V (VA  " V)'  (1) 

P - pA  = pA  (U  - uA)  (u  - uA) , (2) 

where  only  the  one-dimensional  strain  configuration  (particle  or  flow 
velocity  parallel  to  the  shock  velocity) , and  all  quantities  are  measured 
with  respect  to  a fixed  laboratory  coordinate  system.  Additionally,  the 
states  A and  B are  assumed  to  be  in  thermodynamic  equilibrium  so  that  the 
only  forces  present  are  those  associated  with  the  equation  of  state  pressures 
PB  (deonted  by  P)  and  PA . Viscosity  and  thermal  conductivity,  while 
possibly  contributing  to  the  structure  of  the  shock  fronts,  are  assumed  to 
be  unimportant  at  the  positions  where  (P,  u)  are  measured.  This  allows 
the  P values  and  their  derivatives  to  be  characterized  as  "Hugoniot"  values. 
U denotes  shock  velocity. 


Substituting  Eq.  (1)  into  Eq.  (2)  gives 


dP 

dV 


P-PA  = 


u - u , 


(VA  - V) 


H 


vA-v 


2 ^ 
dV 


u - u , 


H + V, 


V 


(3) 


du 


dV  1 H 


a 


(4) 


*When  the  subscript  H indicating  Hugoniot  appears  it  will  indicate  a quantity 
characterizing  a state  in  stable  thermodynamic  equilibrium,  and  that  the  state 
is  achieved  via  the  shock  jump  conditions  (ref.  1)  . 
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In  the  units  of  Eq.  (1)  specific  volume  is  written  as  V = p"1  . Thus,  upon 
rewriting  Eq.  (1) 


u - u. 


U - uA  = 


pa(va-v)  ' 


U = 


u - u, 


V - V/ 


va+ua 


U ~ - VA  I 5V  1 H J + UA 


(5) 

(6) 

(7) 


Equation  (7)  and  the  right  hand  side  of  Equation  (4)  assume  that  only 
small  excursions  from  the  state  A into  the  state  B are  being  observed . 
Thus,  within  the  state  B,  for  an  elastic  constant  M0, 


dP  = - M 


B 


dV 


(8) 


Substituting  Eq.  (8)  intoEq.  (7)  yields 


dP  . _ ^B 

du  1 H u “ UA 


(9) 


dP 

If  the  jump  condition  of  Eq . (2)  holds,  then  — | ^ < 0 for  small 

excursions  from  the  state  A . Further,  from  Eq.  (2),  if  u > uA  and 

dP 

P > PA,  then  U > uA  so  that  — | H > 0. 


If,  however,  large  excursions  from  a known  state  A to  a new  state  B 
are  of  interest,  then  the  ratio  (u  - uA)/(V  - VA)  can  not  be  replaced  by 

first-order  derivatives,  and  more  physics  becomes  necessary  in  order  to 

determine  the  sign  of  | ^ at  the  new  state. 


In  addition  to  being  able  to  make  a similar  statement  for  the  sign  of 
HP 

— I , it  should  be  noted  that  even  small  excursions  to  a new  state 
du  1 H 

allow  for  a sign  change  if  the  new  state  is  characterized  by  (u  < uA  , 
P > PA) , or  (u  > uA,  P < PA. 
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Among  other  considerations,  the  sign  of  I H is  important,  because 
it  bears  directly  upon  the  possibility  of  amplification. 

AMPLIFICATION 

A double  wave  structure  with  an  infinitesimal  amplitude  acoustic  dis- 
turbance overtaking  (and  reflecting  from)  the  second  shock  is  illustrated 
schematically  in  figure  3. 


PB3*  PB3 


pB2>  pb2 


PB1 * PB1 


Figure  3.  Reflection  from  a shock  front. 

Shock  A moving  into  previously  undisturbed  material  (subscript  zero) 
is  followed  by  shock  B which  places  the  medium  into  the  Hugoniot  state 
(p  P ).  A small  amplitude  (acoustic)  compressive  disturbance  orig- 

inating  somewhere  within  (Pg-j/  P^)  causes  a state  change  to  (p^'  PB2^ 

and  then  reflects  from  the  shock  front  B . The  reflected  wave  causes  a 
state  change  from  (PB2'  PB2)  to  ^PB3'  PB3^  ’ The  states  62  and  83  are 
only  infinitesimally  removed  from  the  state  B1 . 

The  following  discussion  of  the  amplification  factor,  (Pg3  - PB2)/ 

(P  -P  ),  is  essentially  due  to  the  work  of  Fowles  (ref.  1).  For 

the  incoming  and  reflected  disturbances,  away  from  the  region  of  re- 
flection. 


r 


(10) 


PB2  " PB1  =1  du  1 s+  ] [UB2  ' UBl]  ' 
PB3  ‘ PB2=(  d^'  s-1  (UB3  “ UB2  ] ' 


where  the  infinitesimal  amplitude  of  the  disturbances  allows  the  assump- 
tion that  the  propagation  processes  are  isentropic  (entropy  changes  are 
third  order  with  respect  to  first-order  pressure  changes  (ref.  5))  . The 
subscript,  s denotes  isentropic,  and  the  (+,  -)  denote  "directionality" 

(see  below)  . 

The  net  change  from  the  state  B1  to  the  state  B3  cannot  be  simply  treated 
as  an  isentropic  process  because  of  a possible  energy  transfer  from  the 
shock  front  to  the  reflecting  disturbance.  It  will  be  argued,  however, 
that  the  state  B3  can  be  achieved  via  a Hugoniot  process  with  respect  to 
the  state  B1 . 

Consider  the  infinitesimal  amplitude  acoustic  disturbance  to  be  in  pulse 
form.  A short  time  after  the  leading  edge  of  the  pulse  is  reflected,  the 
situation  is  as  shown  in  figure  4. 


(CB)si* 83 


Figure  4.  Acoustic  pulse  reflected  from  the  shock  front  B . 

is  the  isentropic  acoustic  velocity  with  respect 
to  the  state  B . 


The  argument  is  simply  that  the  shaded  boundaries  represent  acoustic 
fronts,  and,  that  if  the  state  B3  is  thermodynamically  stable,  it  must  be 
achievable  via  a dynamic  (P,  V,  T)  process.  Thus, 


dP 


du  1 H 


(12) 


Using  equations  (10)  through  (12),  with  the  subscript  B understood, 
forms 

(u3  - u.j ) - (u3  - u2)  - (u2  - Uj)  =0  (13) 


which  results  in 


For  infinitesimal  amplitude  disturbances  propagating  with  respect  to  the 
preexisting  B state,  equation  (2)  becomes 


dp  + = PB  (U+ - ug  ) du  = + PgCg  du. 


(15) 


where  U±  is  the  velocity  of  the  disturbance  front  in  the  laboratory  refer- 
ence frame.  If  Cg  is  the  disturbance  (acoustic)  velocity  with  respect  to  a 

stationary  state  B,  the  U+  = ±Cg  + Ug,  with  Cg  > 0 and  the  positive  sign 

denoting  propagation  in  the  direction  and  sense  of  Ug.  Thus,  upon  com- 
bining equations  (14)  and  (15)  with 


yields 


1 PBCB  dpi H 
1 + PBCB  ~dP  I H 


(16) 


(17) 


It  is  understood  that  the  B1  state  (B  state  for  short)  properties  are  to  be 
used  in  evaluating  the  right  hand  side  of  equation  (17) . 
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If,  in  equation  (17) , ^ | > 0,  as  is  almost  universally  the  case, 

then  { (P3  - P2)/ (P2  ~ P^)}  < 1 and  amplification  upon  reflection  does 

not  occur*.  On  the  other  hand,  if  | < 0,  then  the  right  hand  side 

of  equation  (17)  is  greater  than  unity,  and  amplification  upon  reflection 
occurs. 


There  is  a question  of  the  full  implication  of  the  gain  predicted  by 
equation  (17) . Whether  or  not  one  reflected  pulse  continues  to  grow  in 
amplitude  depends  upon  the  details  of  a further  reflection  from  an  imagined 
boundary  far  to  the  left  of  the  shock  B in  figure  4.  If,  however,  a great 
number  of  pulses  are  generated  within  the  bulk  of  the  B1  state,  then  the 
question  may  become  not  one  of  single  -pass  amplification  but  one  of  a 
net  increase  in  pressure  behind  the  shock  front  B.  The  net  increase  in 
pressure  criteria  would  be  the  more  interesting  requirement  if  the  sub- 
ject were  explosives  modeled  with  pulse  generation  (via  local  energy 
release)  probability  proportional  to  the  pressure. 


The  net  increase  in  the  pressure  requirement  translates  to 
{ (P3  - P2) / (P2  " Pf ) } > -1/  which  from  equation  (17)  is  equivalent 

to**. 


PBCB 


du 

dP 


-iHi  >-i. 


(18a) 


clearly  equation  (18a)  overlaps  the  single-pass  amplification  require- 
ment 


du 

dP 


-'h  <0 


(18b) 


*The  reflected  pulse,  which  eventually  propagates  into  the  state  B1  with 
relative  velocity  Cg_  (fig.  4) , has  an  amplitude  relative  to  the  B1  state 
given  by  (P3  - P^  . 

**From  Eq.  (17)  we  see  that  { (P3  - P2)/(P2  ~ Pi) } > -1  can  be  written  as 
{ (1  - X) /(I  + X) } > -1  with  x =pDcD  —I,,.  This  last  inequality  is 

D D Or  H 

satisfied  by 

(a)  1 - X > -1  - X when  1 + X > 0.  X > - 1 . 

(b)  1 - X < -1  - X when  1 + X < 0.  .‘.No  solution. 


It  is  of  interest  at  this  point  to  anticipate  some  of  the  following  dis- 
cussion in  this  report  and  briefly  inquire  into  the  material  properties 
consistent  with  equations  (18  a and  b) . We  write  the  shock  velocity  - 
particle  velocity  relationship  (ref.  6 and  7) 

U = a + pu 

And,  in  a form  consistent  with  figure  3 

(UB  ' UA)  =<X  + P (UB  " UA}-  (19) 

In  equation  (19),  when  u = uA,  the  shock  velocity  UB  corresponds  to 

an  infinitesimal  disturbance  propagating  into  the  precompressed  state  A. 
Thus,  a = cA  (the  local  sound  speed  in  state  A) . Substituting  Eq.  (19) 
into  Eq.  (2)  yields 


PH  = PACA  (U"UA>  (U'UA>’  +PA  ' 


(20) 


" £lu<0 -Paca  + 2pa  P (u-ua)<0 


du  1 H 
c, 


A A 


p <-2(u-ua) 


(21) 

(22/ 


Thus,  for  u > u , p must  be  less  than  0 in  order  to  satisfy  Eq.  (18  b). 

n 

Huang  (ref.  8)  shows,  subject  to  the  restrictions  of  the  Slater  (ref.  9) 
model. 


{± 
{ ap 


- v 


( lv  1 T ) 1 1 T - » * T • 


(23) 


where  the  subscript  T denotes  the  isothermal  derivative,  that* 

YA^  1 (24) 

2 

Thus,  negative  p implies  a negative  Cruneisen  parameter,  rA,  for  the 
state  A . 


P = -T-  + J 


*There  are  actually  a number  of  different  forms  in  the  literature  for  the 
relationship  between  p and  Ya  (see  references  6 and  7 for  examples 
which  differ  from  Eq.  (24)) . Each  of  those  forms,  however,  leads  to  the 
conclusion  that  negative  p implies  negative  YA. 
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The  Grijneisen  parameter,  y.  is  defined  by 


(25) 


where  E is  specific  internal  energy  and  has  the  equivalent  thermo- 
dynamic definition 


Y = XVJ/Cp  . 


(26) 


X is  the  thermal  expansion  coefficient,  S the  entropy,  Cp  the  specific  heat  at 
constant  pressure,  and  J the  isentropic  bulk  modulus. 


J = - V 


(27) 


J is  related  to  the  isothermal  modulus  by  (ref.  8) 


0V  I s (i  + YX  T)  av  I J • (28) 

dP 

Because  — | ^ < 0 is  required  (ref.  11)  for  thermodynamic  stability 

(with  respect  to  volume  fluctuations) , we  have  J > 0.  Thus,  from  Eq. 

(26) , a negative  y is  associated  with  a negative  thermal  expansion  co- 
efficient. 

Thermal  expansion  measurements  have  found  very  few  materials 
which  exhibit  negative  x-  The  most  common  examples  of  this  unusual 
property  are  silicon  (ref.  12  and  13)  and  indium  antimonide  (InSb) 

(ref.  12)  at  approximately  20°K  and  atmospheric  pressure  and  fused 
(vitreous)  silica  (ref.  12,  14,  and  15)  (Si02)  at  perhaps  room  temperature. 

Silicon  and  InSb  are  not  practical  shock  amplification  candidates,  as 
their  x values  are  only  known  to  be  negative  at  the  inconvenient  temp- 
erature of  20°K  (at  atmospheric  pressure)  . Similarly,  fused  Si02  has, 
at  best,  a small  negative  x at  room  temperature.  What  is  more  important 
than  known  candidates  for  practical  shock  amplification  is  the  simple 
existence  of  these  materials,  and  the  physics  relating  to  that  existence. 

The  physics  picture  (ref.  10)  associated  with  a negative  x involves 
a transverse  phonon  (acoustic  or  optic)  of  relatively  low  frequency.  By 
analogy  with  a stretched  violin  string  (ref.  10)  compression  further  lowers 
the  frequency  (by  reducing  the  violin  string  tension)  . The  single  mode 
Gruneisen  parameter  is  defined  by 
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r-jTsc.--  ^ 


9 In  “i 


V 


(29) 


9 In  V 


oo. 

i 


a (0- 

9 1^ 

9 V 


Clearly,  decreasing  (Oj  in  compression  (decreasing  V)  gives  negative  y. 

As  a general  rule  (ref.  16),  phase  transitions  are  associated  with 
“i  ->•  0 as  the  transition  is  approached . The  thought  then  is  to  utilize 
the  material  properties  associated  with  an  incipient  phase  transforma- 
tion in  order  to  achieve  shock  amplification.  While  the  idea  of  so  utilizing 
the  phase  transformation  processes  was  discussed  previously  (ref.  17 
and  18),  the  approach  as  expressed  in  the  combination  of  Eqs.  (17)  and 
(29)  represents  a new  approach  to  the  physics  of  amplification. 

A negative  mode,  Yj«  of  Ecl-  (29)  does  not  guaranty  that  the  macro- 
scopic y of  Eq.  (26)  is  negative.  The  relationship  between  y and  Yj 
is  complicated;  it  depends  upon  the  detailed  phonon  (mode)  spectrum 
and  the  degree  to  which  each  mode  is  occupied  (i  ,e.  the  degree  to  which 
the  various  lattice  vibrational  states  are  filled)  . Thus,  the  possibility 
of  a negative  Yj  associated  with  phase  transitions  is  a hint  on  where  to 
look  for  candidate  amplification  materials,  rather  than  a prescription 
for  choosing  a specific  material.  It  is  possible  for  Yj  to  t>e  negative  and 
at  the  same  time  to  have  the  macroscopic  y positive.  In  a later  section  of 
this  report  the  question  of  candidate  selection  will  be  treated  in  more 
detail . 


SHOCK  STABILITY  LIMITS 

Analytical  studies  (ref.  19)  have  indicated  that  a shock  disturbance 
which  is  outside  the  limits  given  by 

-1  <j*  (t£ih)  <1+2i^  • (30) 

is  subject  to  an  exponential  temporal  growth*  of  behind  the  shock 
perturbations  of  the  hydrodynamic  parameters. 


*The  analysis  assumes  an  irrotational  scalar  fluid  with  the  growth 
restricted  to  the  region  immediately  behind  the  shock  front. 
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In  Eq.  (30) , M is  the  local  Mach  number  defined  by 

M=  , (31) 

where  c is  the  local  infinitesimal  disturbance  sound  speed  and  j is  the 
mass  flux  relative  to  a coordinate  system  attached  to  the  shock  front*. 

It  is  also  true,  as  shown  by  Fowles  (ref.  1)  that  those  same  limits  of 
Eq.  (30)  correspond  to  the  existence  of  a shock  instability  leading  to 
the  possibility  of  a two  wave  structure.  That  dual  occurrence  of 
exponential  growth  and  shock  instability  represents  an  interesting 
observation  in  that  it  begins  to  connect  a double  wave  structure  with 
the  individual  mode  physics**  ultimately  responsible  for  its  existence. 

Equation  (18)  shows  that  ^ | ^ < 0 is  sufficient  for  either  single 

pass  amplification  or  bulk  pressure  growth.  Thus,  a relationship  be- 
tween du  . and  dV  , is  desired  so  that  the  limits  of  Eq.  (30)  can  be 
dP 1 H dP  1 H 
connected  with  Eq.  (18) . 

From  Equations  (1)  and  (2) 

(u  - uA)2  = (P  - PA)  (VA  - V),  (32) 

/.  2 (u  - uA')du  = (VA  - V)  dP  - (P  - PA)  dV.  (33) 

But,  Equations  (1)  and  (2)  can  also  be  combined  to  give 

p-pa  = '’a'  '%-V  (va-v) 


so  that  jBA  becomes  j^A 


2 _ P " PA 


*For  the  unperturbed  state  B of  figure  3,  and  j become 


MB  = 


'b  - UB 


jBA  PA  (UB  UA} 


**Swan  and  Fowles  (ref.  19)  treat  a linearized  system  so  that  their  results 
may  be  thought  of  as  pertaining  to  a single  mode  of  Fourier  analysis. 


where  it  is  understood  that  (P,  V,  u)  = (P  , V , uR)  . Upon  substituting 
Eq.  (35)  into  Eq.  (32)  B B 

,U-UA>'=1BA  (VA-V)’-  <36) 

Thus,  Eq.  (33)  becomes 

-2jBA(VA"V)  du=  (VA'V)  dP~  (P'PA}  dV'  (37) 


jBA  ( dP  ' h)  -2  1 ' jE 


J dV  | 

BA  dP 


Eq.  (38)  gives  the  desired  connection  between  ^ | and 

Or  H 

dV  | 

dP1  H ’ 

For  a shock  moving  to  the  right,  j . > 0.  Thus,  the  comparison 

BA 

between  Eqs.  (38)  and  (30)  leads  to  the  matrix  shown  in  table  1 . 


Table  1 


Correspondence  between  limit  violations*  for 

du,  . . dV  | 
dpi  h and  dP  I H 


Sign  in  .?  dV  . 

Eq.  (38).  J dP  lH  > 1 + 2 M 


i2  dV,  , 

J dpi  H < " 1 


H<-  J<° 


du  | . -i 

dP  I H > J 


du  | M n 
dP  I H J > 0 


dP  1 H 


.1 

< J < 0 


J and  M again  denote  jRfl  and  MR, respectively, 


Thus,  the  minus  sign  choice  in  Eq.  (38)  presupposes  the  physics 
necessary  for  amplification  as  discussed  in  Eq.  (18),  specifically  in 
Eq.  (18  b)  . Thus,  the  minus  sign  solutions  in  Eq.  (38)  will  be  neglected.* 


The  (P,  u)  and  (P,  V)  plane  possibilities  corresponding  to  the  plus 
sign  choice  are  illustrated  in  figures  5 and  6. 


+ sign  solutions  \ 

B:  Normal  Compression  \ II 

B':  Normal  Rarefaction  . 


J " JBA 


Figure  5.  (P,V)  plane  possibilities  corresponding  to  limit  violations. 

The  chords  (of  slope  -j)  connecting  the  states  A,  B,  and  B’ 
are  drawn  with  double  solid  lines.  The  solid  lines  through 
B,  B‘  states  represent  limit  violating  derivatives  (slopes) 
from  Eq.  (30) . Lines  I and  II  correspond  to  lines  I and  II, 
respectively,  in  figure  6. 

*By  concentrating  on  the  plus  sign  solutions  of  Eq.  (38)  our  investigation 
will  be  restricted  to  positive  chord  B states  (i.e.,  Au  . ) which  con- 

An  ^ " 


currently  admit  negative  local  du  ■ derivatives  (and  thus  result  in 
amplification) . dP  * H 


-j/M 

I 


+ sign  solutions 
B:  Normal  Compression 
B":  Normal  Rarefaction 


Figure  6.  (P,  u)  plane  possibilities  corresponding  to  limit  violations. 

The  solid  lines  through  the  B,  B'  states  represent  limit 
violating  derivatives  from  Eqs.  (30)  and  (38) . Lines  I and 
II  correspond  to  lines  I and  II,  respectively,  in  figure  5. 

While  the  slope  I states  obviously  satisfy  the  ^ | <0  requirement  of 

OP  H 

Eq.  (18  b)  for  single-pass  amplification,  the  situation  with  respect  to  the 
gross  pressure  generation  requirement  of  Eq.  (18  a)  necessitates  analysis 


In  order  for  line  I of  figure  6 to  fit  within  the  inequality  of  Eq.  (18  a) 
it  is  necessary  that  (jAB/pBcB^B)  <1.  But*, 


*While  Eq.  (39)  predicts  U^/pgCgMg)  < 0 (and  thus  also  < 1)  for  a 

shock  moving  to  the  left,  it  is  not  obvious  that  such  behavior  is  to  be 
expected.  The  change  in  sign  for  a wave  moving  to  the  left  is  a con- 
sequence of  the  absolute  value  operator  contained  in  the  definition  of 
Mg • That  absolute  value  appears  in  the  derivation  of  Eq.  (30)  when 

(see  Eq.(38)  of  reference  19  ) an  equation  of  the  form  x < 1 is  replaced 
by  | x ei0  | < 1 and  consequently  may  be  forcing  a directionality  re- 
lated sign  difficulty. 
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JAB 


pbcbmb 


pA(UB 


pbcbI 


V 


UB  UB 


'B 


PB  (UB  ’ UB) 


PBCB 


UB  UB 
CB 


JAB 


pBCBnB 


= 1 


Thus,  violating  the  upper  limit  of  Eq.  (30)  does  not  satisfy  the  inequality 
of  Eq.  (18  a) . Line  II  of  figure  6,  corresponding  to  the  lower  limit  of  Eq. 
(30),  automatically  satisfies  Eq.  (18  a),  as  line  II  is  characterized  by 

-I 

du1  H 


^1  0. 


In  summarizing  this  section,  note  that  Eq . (17),  Hugoniot  require- 
ments for  single-pass  amplification,  or  for  simple  pressure  increase,  are 
each  separately  satisfied  by  one  of  the  limits  of  Eq.  (30)  . Thus,  there  is 
agreement  between  the  detailed  analytic  derivation  (ref.  19)  leading  to 
Eq.  (30)  and  the  less  detailed,  but  conceputally  satifying  picture  associated 
with  Eq.  (18  a and  b)  . 


THERMODYNAMIC  STABILITY 


The  single-pass  amplification  requirement,  dP 


du  1 H 


< 0 


, corresponds 


to  dP 


I..  > 0. 


But,  9P 


9P  . is  required  (ref.  11)  for  thermodynamic 

9V  1 T < ° 


dV  1 H 

stability  with  respect  to  volume  fluctuations.  It  is,  therefore,  necessary 
to  investigate  the  relationship  between  3P  , and  dP, 

3V  I T dV>  H 


=r*.  I . . > 0 


hold  simultaneously,  the 


Should  both  3P  . . . and  dP 

3V  1 T dV  1 H 

implication  would  be  the  existence  of  a stable  shock  state  B which  is 
capable  of  amplifying  behind  the  shock  signals  WHICH  COULD  ORIGINATE 
VIA  VOLUME  FLUCTUATIONS.  Thus,  we  have  a contradiction  and,  in 
turn,  the  implication  that  such  a state  B (in  figs.  1 . 2,  and  3)  could  not 
exist  (in  the  steady  state  sense) , and  that  a new  stable  state  C,  character- 


and  dP 


— < 0 


would  be  the  final  result.  Figure  7 


ized  by  aP  • 

aV  1 T dV  1 H 

illustrates  the  expected  possibilities  in  the  (P,  V)  plane. 


18 


> ■ r 





Figure  7.  Possible  final  stable  states  C2  or  C2  achievable  via 
jjP  . with  respect  to  the  unstable  state  B. 

dV  1 H > ° 

If,  as  illustrated  by  state  C2  in  the  exaggerated  figure  7,  the  cord 
C2  A is  steeper  than  chord  AV0,  then  the  resulting  constant  velocity 
piston  induced  shock  structure  is  as  shown  in  figure  8.  The  shock 
structure  for  the  exaggerated  Cx  state  is  illustrated  in  figure  9.  Both 
C states  are  characterized  by  dP  , < . 

dV  1 H 


C2 


distance 


Figure  8.  Constant  velocity  piston  wave  structure  corresponding  to 
the  exaggerated  state  C2  . 


19 


A 


distance 


Figure  9.  Constant  velocity  piston  wave  structure  corresponding  to 
the  exaggerated  state  Cx  . 

We  will  now  investigate,  in  algebraic  detail,  the  compatibility 

between  3P  . < 0 and  I > 0 ' 

9V  1 T dV  H 

Let  P = P (V,  T)  . Thus, 


dP  , _ 9P  | . £T  | — I 

dV  H 3V  ^ T 9T  V dV  H 

For  a quasi-static  process  within  the  state  B 
TdS  = dE  + PdV . 


(40) 


(41) 


Rewriting  Eq . (41)  as 

) = (dE  | ) + P CdV  | ) « 

p P H 


T Hip  (dT 


where  S,  E.  and  V are  specific  (i.e..  per  unit  mass)  quantities, 
immediately  leads  to 

_ 51  i = — I + p. 

cp  sv  >p  av1  p 

Now,  let  E = E (P,  V).  Then, 


dV  1 H 


2§| 

9V1  P 


|3E,  dP | . 

9P  1 V dV  1 H 


(42) 


(43) 
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Combining  Eqs.  (42),  (43)  and  (25),  yields 


fc  ( - 

1 1 -p] 

V 

i 

( - 1 1 

Lcp  1 av 

1 pi  PJ 

T 

y 

1 dV  ' H ) 

V) , then 

51  1 + 

51, 

av 1 p 

ap1  v 

1 dV  1 

H ‘ 

Because  dE  , is  a known  function  of  the  shock  structure  (see 

dV  H dT 
below),  Eqs.  (44)  and  (45)  can  be  utilized  to  eliminate  ^ | H from 

Eq.  (40).  H is  then  given  in  terms  of  ~|  T and  other  thermo 

dynamic  derivatives  and  the  observational  parameters  associated  with 
the  shock  state.  Combining  Eqs.  (44)  and  (45)  gives 


dT  | J_  dE  , 

dV  ' H c dV  1 H 


51 1 + -¥-  5E  | 

3P 1 V cnY  / dV  H c 


Upon  substituting  Eq.  (46)  into  Eq.  (40), 

V 3P  , \ dP  , 

| 2 + CpY  3T  V I dV  H 


aP,  + _L  — | — i + -L. 

3V1  T cp  1 3T  1 VI  dV  1 H cp 


3T  1 V 


The  Hugoniot  energy  relation,  for  the  state  B propagating  into  the 
state  A,  has  the  form  (ref.  1,  7) 


E - EA  = i (p  + FV  (Va-v)- 


••  - I 

dV  1 H 


dP. 

2 2 dV  1 H 


Substituting  Eq.  (49)  into  Eq.  (47)  yields,  after  minor  rearrangement. 


Eq.  (50)  is  the  desired  expression  for  the  relationship  between  ~ I .. 
ar>d  gy  I j • With  the  exception  of  (^A  and  VA) , all  quantities  on  the 
right  hand  side  of  Eq.  (50)  are  to  be  evaluated  at  the  state  B . 

A sign  analysis  of  Eq.  (50)  is  performed  next,  assuming  tt.  I _ < 0 

oV  T 

for  thermodynamic  stability.  From  reference  20 


-I 

3T  1 V 


— I 

av 1 p 


-I 

av 1 t ' 


and  Eqs.  (26),  (27),  and  (28),  it  follows  that 


(51) 


9P, 

P l 3T  1 V 


XVJ 

Y 


,-1 

av 1 

_ VQ  + yxT ) 


3T 

5V 


-X 


-I 

av 1 t 


(52) 

(53) 


Substituting  Eq.  (53)  into  Eq.  (50)  for  AP  and  AV,  each  being  small 
and  of  arbitrary  sign,  yields 


. dP  , _ 3P  , 

'■  dV  I H ~ av  I T 


L±JOL 
. 2 + yxt 


(54) 


dP  ftP 

Thus,  for  small  AP  and  AV»  I H < 0 if  — | < 0 and  | yxT|  < 1 

REGARDLESS  of  the  sign  of  y.  Thus,  thermodynamic  stability  requires 
dP 

— | H < 0 regardless  of  the  sign  of  y for  small  pressure  and  volume 
excursions  from  the  stable  state  A. 
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It  is  possible  to  take  the  sign  analysis  a step  further.  By  imposing 
J > 0 as  a requirement*,  and  employing  the  extreme  right  hand  side  of 
Eq.  (52)  in  Eq.  (50),  the  numerator  within  the  square  brackets  of  Eq. 
(50)  will  be  negative  if  y < 0,  P > P a,  and  9P  , < . Since  the 

9V  1 T 

denomenator  within  the  square  brackets  of  Eq . (50)  can  be  written 


lap.  | *x 

—i 

VA  - V (1  + 2y  ) 

Y 

\ sv'ti 

2 

dP 

certainly  ^ | ^ < 0 will  hold  if 
inequality  implies 


VA  - v (1  + \ ) 


> 0. 


This  last 


1 ^ VA  - V 

y 2 V 

But,  we  have  already  chosen  y<0  so  that  Eq.  (55)  becomes 


(55) 


Irl  > v^T' 

which  requires  y to  be  large  and  negative. 


(56) 


The  above  paragraph  shows  that  even  for  a large  negative  y, 

— I < 0 if  |r,  I t < 0 holds.  Thus,  the  requirement  of  thermodynamic 
dv  1 H 3V  T 

stability  (with  respect  to  volume  fluctuations)  appears  to  restrict  dP  . 

d V 1 H 

to  values  such  that  the  upper  inequality  of  Eq.  (30)  is  not  violated.  The 
discussion  associated  with  Eq.  (14)  showed  that  dP  . < n is  associated 

du  1 H 

with  y<  0.  The  implication  is  clear;  dP  , . js  inconsistent  with  a 

dV  1 H 

state  of  thermodynamic  equilibrium.  Consequently,  a state  characterized 
by  dP  - is  expected  to  be  metastable. 

dV  ' H ° 


9P 

*ln  order  that  a real  local  adiabatic  sound  velocity  given  by  c2  = * V2  — 

exist. 
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What  does  all  of  this  mean?  If  the  sign  analyses  were  correct,  it 
would  imply  that  y<  0 materials  do  not  exist  in  equilibrium,  and  that  a 
medium  capable  of  supporting  single  pass  amplification  must  be  formed 
into  a metastable  state  (with  some  finite  lifetime)  before  that  single-pass 
pressure  amplification  could  be  expected  to  occur.  Correspondingly,  if 
one  were  to  identify  the  dP  | > Q limit  violation  of  Eq.  (30)  with  a 

detonation  (as  does  Fowles,  reference  1),  then  the  implication  would  be 
that  there  is  little  utility  in  studying  the  physics  of  explosives  which  are 
in  an  equilibrium  state  A.  The  physical  properties  of  the  y<  0 non- 
equilibrium state  B must  be  so  radically  removed  from  the  state  A that  in 
effect  one  would  be  dealing  with  two  radically  different  (isomorphic) 
materials . 


However,  it  is  known  that  states  of  macroscopic  y<  0 (e.g..  Si  and 
InSb  at  low  temperatures)  exist  in  apparent  thermodynamic  equilibrium'. 
This  is  a puzzling  paradox. 

MATERIAL  SELECTION 


In  this  section,  the  possibility  of  finding  non-energetic  materials 
with  y<  0 in  the  vicinity  of  room  temperature  is  investigated.  The 
association  between  an  individual  mode  (Yj<  0)  and  phase  transitions 
was  already  briefly  discussed . Now  the  y-phase  transition  relationship 
will  be  considered  from  a macroscopic  point  of  view. 

Two  classes  of  phase  transitions  can  be  identified  as  depending 
upon  V vs  T on  either  side  of  the  transition.  The  situations  are  illustrated 
in  figures  10  and  11,  along  with  the  corresponding  thermal  expansion  co- 
efficients x- 


CLASS  I 


/ ; 

i 

| I 

TC'  T 

Figure  10.  Phase  transition  resulting  in  x > 0.  Tc  denotes  the  transition 
temperature.  The  thermal  expansion  coefficient  is  displayed 
by  the  dashed  curve. 
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Figure  1 1 . Phase  transition  resulting  in  x < 0-  Tc  denotes  the  transition 
temperature.  The  thermal  expansion  coefficient  is  displayed 
by  the  dashed  curve. 

T ixN i i _x , depending  upon  heat  treatment  and  composition,  is  capable 

of  exhibiting  either  Class  I or  Class  II  behavior  in  the  process  of  under- 
going the  TiNi  (II  ■*-*■  III)  phase  transition  (ref.  21).  Class  I behavior  is 
exhibited  (ref.  22)  by  approximately  stoichiometric  TiNi  and  Class  II 
behavior  by  (ref.  22)  TiNi  compounds  containing  approximately  55  percent 
nickel . 

The  important  observation  about  the  possibly  negative  macroscopic  y 
associated  with  the  Class  II  behavior  illustrated  in  figure  11  is  that  y < 0 
is  to  be  found  only  in  a narrow  temperature  band  surrounding  Tc.  That 
observation  supports  the  tentative  conclusions  (see  the  discussion  asso- 
ciated with  the  inequality  of  Eq.  ( 56) ) in  that  a y < 0 state  implies  meta- 
stability rather  than  a long-lived  state  in  thermodynamic  equilibrium. 


Thus,  a Class  II  material  could  not  be  used  for  continuous  single- 


pass amplification  (in  the  sense  of  dP 


lu<0 


) , because  volume  fluctua- 


du  1 H 

tions  behind  the  shock  front  B would  tend  to  trigger  the  phase  transition 
to  a y > 0 state  on  one  side  of  Tc . 
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CONCLUDING  DISCUSSION 


Figure  12  illustrates  an  exaggerated  view  of  a desired  result  of 
single-pass  pressure  amplification. 


before  interaction 


J-L 


immediately  after  interaction 

r : 

! i B + C' 


A 

i 


A 

1 


Figure  12.  A finite  amplitude  pulse  C introduced  into  the  state  B and 
overtakes  the  shock  front  B.  Upon  total  reflection  from  the 
front  B its  amplitude  is  increased,  (C) . The  width  of  C' 
is  shown  as  less  than  the  width  of  C in  order  to  conserve 
strain  energy. 

Where  it  is  possible  to  extract  the  pulse  C1  of  figure  12  and  introduce 
it  into  an  explosive,  one  would  have  a more  efficient  means  of  explosive 
initiation  in  the  sense  of  the  P2t  criteria  (ref.  2),  as  discussed  earlier . 


Unfortunately,  previous  arguments  appear  to  require  that  an  inert 
amplifying  state  B (macroscopic  y<  0)  be  metastable  and  possibly  associated 
with  an  occurring  phase  transition  (microscopic,  or  mode,  Yi<  • Because 
the  state  B is  metastable,  it  cannot  be  maintained  while  waiting  for  the 
introduction  of  the  pulse  C.  Thus,  practical  single-pass  pressure  amp- 
lification does  not  seem  to  be  possible. 


The  above  argument  on  inert  state  inaccessibility  due  to  metastabil  ity 
does  not  hold  for  an  energetic  (e.g.  explosive)  medium*.  Indeed,  it  is 
known  (ref.  23)  that  an  ongoing  energetic  chemical  reaction  can  serve  as 
an  amplifying  medium  for  acoustic  signals.  Such  an  effect i can  ‘ be  though 
of  as  being  allowed  basically  because  thermal  energy  (release)  per  unit 


*Equation  (48)  must  be  modified  to  take  energy  release  into  account  in 
treating  energetic  materials. 
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volume  is  dimensionally  equivalent  to  pressure  (increase) . Thus,  if 
propagation  conditions  (e.g.  Chapman  - Jouguet  conditions  (ref.  24)) 
are  fulfilled,  an  induced  metastability  leading  to  energy  release  (pres- 
sure increase)  finds  itself  in  phase  with  an  amplified  propagating 
disturbance. 
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